An efficient Bayesian inference approach to inverse problems based on adaptive sparse grid collocation method

نویسندگان

  • Xiang Ma
  • Nicholas Zabaras
  • Frank H. T. Rhodes Hall
چکیده

A new approach for modeling inverse problems using a Bayesian inference method is introduced. The Bayesian approach considers the unknown parameters as random variables and seeks the probabilistic distribution of the unknowns. By introducing the concept of the stochastic prior state space to the Bayesian formulation, we reformulate the deterministic forward problem as a stochastic one. The adaptive hierarchical sparse grid collocation method (ASGC) is used for constructing an interpolant of the solution to this stochastic forward model in this prior space which is large enough to capture all the variability/uncertainty in the posterior distribution of the unknown parameters. This solution can be considered as a stochastic function of the random unknowns in the stochastic space. This function serves as a stochastic surrogate model for the likelihood calculation. Hierarchical Bayesian formulation is used to derive the posterior probability density function (PPDF). The spatial model is represented as a convolution of a smooth kernel and Markov random field to achieve dimension reduction for facilitating the inference algorithm. Then the state space of the PPDF is explored using Markov chain Monte Carlo (MCMC) algorithms in order to obtain statistics of the unknowns. Instead of repeated evaluations of the deterministic forward model, the likelihood calculation is reduced by directly sampling the approximate stochastic solution obtained through the ASGC method. The technique is assessed on two non-linear inverse problems: source inversion and permeability estimation in flow through porous media.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

An efficient Bayesian inference approach to inverse problems based on an adaptive sparse grid collocation method

A new approach to modeling inverse problems using a Bayesian inference method is introduced. The Bayesian approach considers the unknown parameters as random variables and seeks the probabilistic distribution of the unknowns. By introducing the concept of the stochastic prior state space to the Bayesian formulation, we reformulate the deterministic forward problem as a stochastic one. The adapt...

متن کامل

A Bayesian approach to multiscale inverse problems using the sequential Monte Carlo method

A new Bayesian computational approach is developed to estimate spatially varying parameters. The sparse grid collocation method is adopted to parameterize the spatial field. Based on a hierarchically structured sparse grid, a multiscale representation of the spatial field is constructed. An adaptive refinement strategy is then used for computing the spatially varying parameter. A sequential Mon...

متن کامل

An adaptive sparse-grid high-order stochastic collocation method for Bayesian inference in groundwater reactive transport modeling

[1] Bayesian analysis has become vital to uncertainty quantification in groundwater modeling, but its application has been hindered by the computational cost associated with numerous model executions required by exploring the posterior probability density function (PPDF) of model parameters. This is particularly the case when the PPDF is estimated using Markov Chain Monte Carlo (MCMC) sampling....

متن کامل

Markov-tree Bayesian Group-sparse Modeling: Efficient Solutions to Large Inverse Problems

In this paper, we propose a new Markov-tree Bayesian modeling of wavelet coefficients. Based on a group-sparse GSM model with 2-layer cascaded Gamma distributions for the variances, the proposed method effectively exploits both intrascale and interscale relationships across wavelet subbands. To determine the posterior distribution, we apply Variational Bayesian inference with a subband adaptive...

متن کامل

A Stochastic Collocation Approach to Bayesian Inference in Inverse Problems

We present an efficient numerical strategy for the Bayesian solution of inverse problems. Stochastic collocation methods, based on generalized polynomial chaos (gPC), are used to construct a polynomial approximation of the forward solution over the support of the prior distribution. This approximation then defines a surrogate posterior probability density that can be evaluated repeatedly at min...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2009